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METHOD AND SYSTEM FOR COORDINATE TRANSFORMATION TO MODEL 
RADIAL FLOW NEAR A SINGULARITY 

RELATED INFORMATION 

This application claims priority under 35 U.S.C. 119(e) to 
provisional patent application No. 60/215,697, filed June 
29 f 2000, entitled "Method and System for Oil Reservoir 
5 Simulation and Modeling," which is hereby fully incorporated 
by reference. 

TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to methods and systems for 
10 modeling physical systems using Finite Element analysis and",, 
more specifically, to methods and systems for modeling 
physical systems near a singularity. Even more particularly, 
this invention relates to a method and system for coordinate 
transformation to model radial flow near a singularity. 

15 

BACKGROUND OF THE INVENTION 

Physical systems can be modeled mathematically to simulate 
their behavior under different conditions. A wide variety 
of means exist to model physical systems, ranging from the 

20 very simplistic to the extremely complicated. One of the 

more complicated means to model physical systems is through 
the use of finite element analysis. As the name implies, 
finite element analysis involves the representation of 
individual, finite elements of a physical system in a 

25 mathematical model and the solution of this model in the 
presence of a predetermined set of boundary conditions. 

In finite element modeling, the region that is to be 
analyzed is broken up into sub-regions called elements. 
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This process of dividing the region into sub-regions may be 
referred to as discretization or mesh generation. The region 
is represented by functions defined over each element. This 
generates a number of local functions that are much simpler 
5 than those which would be required to represent the entire 
region. The next step is to analyze the response for each 
element. This is accomplished by building a matrix that 
defines the interaction properties of the various elements 
within the region and a vector that defines the forces 

10 acting on each element in the structure. Once all the 

element matrices and vectors have been created, they are 
combined into a structure matrix equation. This equation 
relates nodal responses for the entire structure to nodal 
forces. After applying boundary conditions, the structure 

15 matrix equation can be solved to obtain unknown nodal 

responses. Intra-element responses can be interpolated from 
nodal values using the functions which were defined over 
each element . 

As indicated above, finite element modeling involves the 
20 creation of a mesh of finite elements. The elements are 
defined by nodes within the problem space. The nodes are 
simply points in space. The lines between the nodes are 
referred to as "edges." The mesh is typically a structured 
mesh. In other words, the mesh is defined in three 
25 dimensions so that the elements within the problem space are 
hexahedrons. For example, they may be cubes or rectangular 
prisms. (Equivalently, in two dimensional problems, the 
elements would be rectangles - see FIGURE 1 for an 
illustration of these elements.) The edges of the 
3 0 hexahedrons are coincident with the edges between the nodes 
of the mesh. In a simple model, the nodes of the mesh may 
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be regularly spaced to define cubic elements. It is not 
necessary, however, for the nodes to be evenly spaced in all 
finite element models. A variety of different sizes and 
shapes of hexahedral elements can be defined within a single 
5 mesh. 

However, mathematical models to simulate the behavior of 
physical systems can be both time consuming and complex to 
create. In particular, finite element analysis problems in 
multiple dimensions (for example, the four dimensions of 

10 space and time) are almost intractable to solve using non- 
computer-aided computational techniques. These types of 
problems are so complex that they must be performed with 
computer assistance. In particular, the difficulty of 
solving the mathematical equations used to model a physical 

15 system increases almost exponentially, in terms of the 

actual computations that must be performed, as the number of 
dimensions that the problem encompasses increases. 

One method of solving the equations involved in mathematical 
modeling of physical systems, such as oil field reservoirs, 

20 involves the use of simulators that can create a model of 

the physical system based on user inputs of both equations, 
and variables within the equations, that describe the 
system. Typically, analysts skilled in both high-level 
mathematics (to generate the complex equations needed) , and 

25 in high-level programming (to generate the necessary code 
from the mathematical formulations) are required. 
Alternatively, teams of specialized mathematicians and 
programmers can be used to analyze a problem, formulate the 
mathematical equations, and write the corresponding code by 

30 hand. One such method and system for creating simulators to 
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model physical systems is disclosed in related U.S. Patent 
Application Serial No. 09/419,717, entitled "Method and 
System for Generating Software Code Using a Symbolic Language 
Translator," filed on October 14, 1999. This application is 
5 hereby incorporated by reference in its entirety. 

However, even the models obtained with these state-of-the- 
art simulators do not always adequately represent the 
physical system being modeled such that accurate solutions 
can be obtained at certain critical points. For example, 
10 prior art methods and systems for solving the equations used 
to model a physical system typically use linear 
approximations and other averaging techniques (that make 
various assumptions that may not be true at certain points 
in the system) . Also, these prior art methods and systems 
15 typically require a large number of nodes (solution points) 
to accurately model a system. Solving the equations that 
represent the system at each of these nodes requires a great 
deal of time and computational power. 

Further, certain physical characteristics of a system being 
modeled do not lend themselves well to modeling in a 
Cartesian coordinate system. In particular, the pressure 
profile in the vicinity of a well tapping into an oil field 
reservoir is highly non-linear in Cartesian coordinates. 
Therefore, using linear elements in the finite element 
analysis may introduce significant interpolation errors. 
This is especially true when a coarse mesh (as known to 
those in the art) , and therefore large linear elements, is 
used in Cartesian space. If the modeling is performed using 
only the X,Y,Z coordinate (Cartesian coordinate) system, 
then if the temporal (time) dimensional extrusion used in 
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the system modeling is irregular (i.e., the system model is 
not uniform) , then the solution to the problem may 
oscillate. 

Modeling of a physical system is a predictive exercise. For 
5 example, when generating a model of the pressure 

distribution inside of an oil field reservoir, the intent is 
to generate a tool that can be used to determine the 
pressure distribution, under a certain set of relevant 
conditions, to be able to predict the effect on the system 
10 when those conditions are changed. The model must therefore 
be accurate within the range of conditions of interest. 

The pressure distribution in an oil field reservoir has been 
found to be non-linear in a Cartesian coordinate system (in 
fact, it is a logarithm), especially in the vicinity of a 

15 well. The pressure is thus typically greater far out from 

the well and then drops logarithmically as the distance from 
the well bore decreases. Modeling such a system accurately, 
using a Cartesian coordinate system, thus requires a large 
number of nodes. A linear approximation is insufficient, 

20 because on the logarithmic pressure curve the values will 
change too greatly between nodes. Using a Cartesian 
coordinate system can thus result in oscillations and 
inaccuracies in the solution. Without an appropriate 
approximation,, predictive behavior is uncertain and the 

25 model fails. 

Standard finite element analysis techniques applied to the 
Cartesian coordinate space are thus inadequate for modeling 
radial flow near a singularity such as the near-well region 
of an oil field reservoir. Boundary conditions representing 
30 fluid flows into, or out of, a well are inaccurate near such 
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a singularity unless a very fine spatial mesh requiring a 
large number of nodes is generated. Such a mesh requires 
nodes spaced at very small increments to be able to 
accurately predict results. The time and computational 
5 resources needed to obtain a solution in this matter also 
increase as the node count increases . 

SUMMARY OF THE INVENTION 

Therefore, a need exists for a coordinate transformation 
10 method and system for modeling radial flow near a 

singularity that can transform the equations used to model a 
physical system from the Cartesian coordinate system to a 
coordinate system in which the solutions to the transformed 
equations are linear for a quantity of interest, such as the 
15 pressure distribution within a reservoir. 

A further need exists for a coordinate transformation method 
and system for modeling radial flow near a singularity that 
can yield at least as accurate solutions as prior art finite 
element analysis techniques applied in a Cartesian 
2 0 coordinate space while requiring a lesser number of nodes to 
model a system. 

Still further, a need exists for a coordinate transformation 
method and system for modeling radial flow near a 
singularity that can provide greater modeling accuracy and 
25 computational efficiency than currently existing finite 

element techniques applied in a Cartesian coordinate space. 

Even further, a need exists for a coordinate transformation 
method and system for modeling radial flow near a 
singularity that can be used together with a non- transformed 
30 Cartesian coordinate model to provide accurate solutions 
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when modeling a physical system near a singularity and in 
the regions spreading outward from the singularity. 

In accordance with the present invention, a coordinate 
transformation method and system for modeling radial flow 
5 near a singularity are provided that substantially eliminate 
or reduce the disadvantages and problems associated with 
currently existing finite element analysis systems and 
methods applied in a Cartesian coordinate space. In 
particular, the present invention provides a method and 

10 system for predicting the behavior of a physical system. 

One embodiment of the method of this invention comprises the 
steps of creating an equation in a first coordinate system 
to model an aspect of the physical system; applying a 
coordinate transformation to the equation to transform the 

15 equation from the first coordinate system into a second 
coordinate system more closely representative of an 
analytical solution to the equation; solving the equation in 
the second coordinate system to obtain a solution ; 
transforming the solution back to the first coordinate 

20 system; creating a second equation in the first coordinate 
system to model a second aspect of the physical system; 
solving the second equation in the first coordinate system 
to obtain a solution to the second equation; and combining 
the mapped solution to the first equation and the solution 

25 to the second equation in the first coordinate system to 

obtain a combined solution. One embodiment of the system of 
this invention comprises a computer- readable medium 
containing a plurality of instructions embodying the above 
method . 
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The present invention provides an important technical 
advantage of a coordinate transformation method and system 
for modeling radial flow near a singularity that can be used 
to transform the equations used to model a physical system 
5 from the Cartesian coordinate system to a coordinate system 
in which the solutions to the transformed equations are 
linear for a quantity of interest, such as the pressure 
distribution within a reservoir . 

Further, the present invention provides an important 
10 technical advantage of a coordinate transformation method 

and system for modeling radial flow near a singularity that 
can yield as accurate solutions as prior art finite element 
analysis techniques applied in a Cartesian coordinate space 
while requiring a lesser number of nodes to model a system. 

Even further, the present invention provides an important 
technical advantage of a coordinate transformation method 
and system for modeling radial flow near a singularity that 
provides greater modeling accuracy and computational 
efficiency than currently existing finite element analysis 
techniques applied in a Cartesian coordinate space. 

Still further, the present invention provides an important 
technical advantage of a coordinate transformation method 
and system for modeling radial flow near a singularity that 
can be used together with a non- trans formed Cartesian 
coordinate model to provide accurate solutions near a 
singularity and in the regions spreading outward from the 
singularity. 
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BRTEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the present invention and 
the advantages thereof may be acquired by referring to the 
following description, taken in conjunction with the 
5 accompanying drawings in which like reference numbers 
indicate like features and wherein: 

FIGURE 1 is an illustration of the rectangular and 
hexahedral elements that comprise 'structured finite element 
meshes in the prior art. 

10 FIGURE 2 shows pressure distribution in the vicinity of a 
producing well, plotted against radial distance from the 
well ; 

FIGURE 3 illustrates pressure distribution in the vicinity 
of a producing well, plotted against the natural logarithm 
15 of the distance from the well; 

FIGURES 4 and 5 illustrate pressure distribution in the 
vicinity of a producing well, computed using the embodiments 
of the radial transformation method of the present 
invention; 

FIGURE 6 illustrates pressure distribution in the vicinity 
of a producing well, plotted as a function of radial 
distance from the well, computed using (i) analytical 
solutions, (ii) Cartesian coordinates, and (iii) an 
embodiment of the radial transformation method of the 
present invention; and 



WO 02/03263 



10 



PCT/US01/20814 



FIGURE 7 is an enlarged graph of the same solutions shown in 
FIGURE 6, showing details of the pressure distribution out 
to a distance of 20 feet. 

5 DETAILED DESCRIPTION QF THE INVENTION 

Preferred embodiments of the present invention are 
illustrated in the FIGURES , like numerals being used to 
refer to like and corresponding parts of the various 
drawings . 

10 The various embodiments of the present invention provide a 
coordinate transformation method and system for modeling 
radial flow near a singularity that is more flexible, more 
computationally efficient, more accurate, and less dependent 
on having a large number of nodes to model the system than 

15 prior art finite element analysis methods and systems. In 

particular, the embodiments of the present invention provide 
the capability for taking an equation representing a system 
to be modeled (or some aspect of the system to be modeled) 
in the Cartesian coordinate system and performing a 

2 0 transformation to represent the equation in some other 

coordinate system in which the solution is in a more usable 
form. 

The embodiments of the present invention can be used to 
perform a mapping from one coordinate space to the other to, 

25 for example, provide a representation in the new coordinate 
space that yields a linear solution. The solution to an 
equation that in the Cartesian coordinate system is non- 
linear, such as the pressure distribution in an oil field 
reservoir, can thus be placed in a map to a coordinate 

30 system in which an accurate predictive model that is linear 
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can be generated. Because the non- linear solution is mapped 
into a coordinate system in which the solution is linear, 
accurate solutions can be obtained using a smaller number of 
data nodes (solution points at which various parameters must 
5 be known to solve the equation (s) modeling a system) . 

The embodiments of the method and system of the present 
invention can use a combination of Cartesian operators and 
non-Cartesian operators to generate accurate and usable 
solutions for equations representing a physical system. 

10 Cartesian operators are used to solve for a solution of the 
equations in the regions of the physical system not near a 
singularity. Non-Cartesian operators are transformed into a 
new coordinate system and solved in the transform space. 
The non-Cartesian solutions to the modeling equations are 

15 then mapped back into the Cartesian coordinate system. 

The physical system being modeled (for example, the region 
in the vicinity of a well bore of an oil field reservoir) 
can thus be represented more accurately using one set of 
operators near a singularity (e.g., the well bore) and a 

20 different set of operators in the regions further away from 
the well bore. In the regions very near the well bore, for 
example, to a radius of about 50 to 100 feet from the well 
bore (singularity) , non-Cartesian operators can be used with 
the radial transformation applied to them. The portions of 

25 the physical system extending beyond this distance can be 
modeled using the Cartesian operators. The distances away 
from a singularity at which the Cartesian and non-Cartesian 
operators are used can be set as required for a particular 
application. 
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The combination of Cartesian and non-Cartesian operators 
yields a good solution throughout the physical system (for 
example, an oil field reservoir) being modeled and can 
provide greater accuracy in a critical area of interest, 
5 e.g., a singularity such as a well bore. As can be seen in 
the example of FIGURE 2, the radial pressure near a well 
bore can change greatly over small distances, but flattens 
out and can be linearly approximated starting several 
hundred feet out from the well bore. In the critical area 

10 of interest near the well bore, applying a coordinate 

transformation in accordance with the embodiments of this 
invention can yield greater accuracy using less nodes and 
less computational resources than with prior art finite 
element analysis techniques. In the regions several hundred 

15 feet out from the well bore, the relationship between 

pressure and distance varies much less rapidly and can be 
more easily approximated using standard techniques. To 
obtain the same accuracy as possible with the embodiments of 
the method and system of the present invention, prior art 

2 0 systems required a much greater number of data nodes in the 

near- well region to provide an accurate model. This is due 
to the rapidly changing pressure profile over short 
distances of this near-well region. 

Under certain circumstances, the radial coordinate 
25 transformation provided by the embodiments of the present 

invention can provide an exact solution. For example, if an 
oil field reservoir being modeled is homogeneous, it may 
only take two nodes to calculate a solution out to any 
radius. Intermediate nodes are not necessary, and the node 

3 0 count required to obtain an accurate solution is greatly 

reduced, resulting in greatly increased computational speed 
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and in a reduction in computational resources required to 
obtain a solution. In such a case, the coordinate 
transformation operators can be used to obtain a solution 
throughout the oil field reservoir model, without the need 
5 for the Cartesian operators in the regions further from the 
well bore (singularity) . 

A lack of resolution in accuracy in the near-well region can 
also arise when using Cartesian discretization in the near- 
well region. A very high mesh resolution, and therefore a 

10 large node count, is necessary to achieve sufficient 

accuracy when compared with a reference analytical solution 
when using the Cartesian operators. The embodiments of the 
method and system for coordinate transformation of the 
present invention can greatly reduce or eliminate both the 

15 lack of resolution in the accuracy in the near-well region 
and the oscillations due to irregular time extrusions. 

The embodiments of the method and system of this invention 
take into account that the flow field around the well in an 
oil field reservoir is essentially radial. The pressure 

20 profile is highly non-linear with respect to distance from 
the well-bore and is poorly approximated by linear elements. 
Radial flow implies that pressure varies linearly against 
the log of radial distance. The embodiments of the present 
invention take advantage of this relationship to 

25 significantly improve the pressure field interpolation. 

This is illustrated in FIGURES 2 and 3, which show a radial 
pressure profile versus distance from a well bore in both a 
Cartesian coordinate system (FIGURE 2) and using the 
transform operators of the embodiments of the present 

30 invention (FIGURE 3) . 
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A basis for the coordinate transformation method of the 
embodiments of the present invention is that steady state 
radial flow from a reservoir towards a singularity (e.g., a 
well bore) is described by the following relationship: 

5 

Q= S^) (P - P w ) <EQH - " 

In Equation 1, W Q" is the flow rate, tt p" is the pressure at 
a reference radius "r" and the subscript w w" denotes the 
10 well bore (singularity) . By rearranging Equation 1, the 
pressure distribution for a given flow rate tt Q" and well 
pressure w p w " is found to be: 



15 



P-P + 



QMr/r w ) 
(EQN. 2) 



The relationship of Equation 2 indicates that a plot of 
pressures versus the natural log of "r" will yield a linear 
relationship. Equations 1 and 2 assume unit conductivity in 

20 the reservoir. FIGURE 2 is a plot of the radial pressure 

profile of a reservoir in Cartesian coordinates. As can be 
seen in FIGURE 2, pressure falls off rapidly near the well 
bore (the first several hundred feet out from the well) . 
Pressure varies much less rapidly the greater the distance 

25 from the well bore and then tends to flatten out. The same 
data plotted on a logarithmic horizontal axis yields a 
straight line relationship, as shown in FIGURE 3. FIGURE 3 
demonstrates that steady state radial flow implies that 
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pressure varies linearly with the log of radial distance 
from a well bore. 

Because the pressure profile for an oil reservoir is highly 
5 non-linear in Cartesian coordinates in the vicinity of the 
well bore (e.g., for distances less than two hundred fifty 
feet in the example shown in FIGURE 2) , using linear 
elements in the finite element analysis may introduce 
significant interpolation errors. This is especially true 

10 when a coarse mesh, and therefore large linear elements, is 
used in Cartesian space. The linear relationship between 
the radial pressure distribution and the log of radial 
distance is therefore a basis for performing the coordinate 
transform. By transforming from the Cartesian (x, y, z) 

15 coordinates into a modified form of radial coordinates (s, 
0, z) where: 

s = ln(r) (EQN. 3) 

20 The embodiments of the method and system of this invention 
can take advantage of the fact that, in the region where 
radial flow is predominant, pressure varies linearly with 
"s" . Therefore, when linear elements are used, and when the 
nodal values are exact, there will be no interpolation 

25 error. 

The following equations and explanation provide more details 
on the embodiments of the coordinate transformation method 
and system of this invention, as applied to the mass 
30 conservation equation for single-phase flow. These 

equations can be used to describe radial fluid flow near a 
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singularity, such as a well bore. The continuity equation 
for single-phase flow is a diffusivity equation of the form 

%-V 2 p = 0 (EQN. 4) 

dt 

The Laplacian operation (V^, in" radial coordinate form, is: 

d 2 p 1 d f dp} 1 d 2 p d 2 p 

2 r ^T + "7T# + -T (EQN. 5) 

Applying the transform s=ln(r), the differentials can be 
written as: 



or ds 



(EQN. 6) 
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The 



continuity equation can be rearranged as: 



f-Exp[-2s] 



yds 2 dS 2 j 



-2 

-f-£. = 0 (EQN. 7) 
6z 2 



10 



In applying the transformation to the one-phase black-oil 
equation, permeability anisotropy is accounted for by 
properly accounting for the cross terms introduced by the 
transform. Given a 3x3 permeability tensor, K, it is 
transformed to radial coordinate form by the following 
expression. 



cos0 sin# 0 
R= -sin# cosd 0 
0 0 1 



The embodiments of the coordinate transformation to model 
15 radial flow near a singularity of this invention can thus be 
used to significantly improve the modeling of one-phase 
fluid flow. 

FIGURES 4 and 5 are graphs illustrating the radial fluid 
flow characteristics around a singularity using the 

20 coordinate transformation model of the embodiments of the 
present invention. The system modeled in this case, as 
shown in FIGURE 4, is a one-quarter segment of a well in the 
center of a circular, reservoir . The reservoir radius is 
1,500 feet, and a relatively coarse mesh is used in the 

25 vicinity of the well (singularity) in the foreground of 

FIGURE 4. The inner radius is equal to .25 feet, the outer 
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radius is 100 feet, with four intermediate computation nodes 
between these two radii. Boundary conditions are that the 
reservoir is initially at 2,000 psi and a constant 
withdrawal is imposed at the well at all times (from start 
5 to finish at ten days) . 

FIGURE 5 is a graph of well pressure versus time, showing 
the analytical solution 40, Cartesian operator solution 50, 
and a transform coordinate solution 60. The plot shown in 
FIGURE 5 comprises forty time intervals, which can help to 

10 improve the overall accuracy of the solution (accuracy can 
improve with the number of time intervals used) . As can be 
seen in FIGURE 5, accumulated computational errors toward 
the end of the computed time show up very clearly in the 
Cartesian operator solution 50, but not in the radial 

15 (transform) coordinate solution 50. The transform 

coordinate solution 60 also tracks much more closely the 
observed analytical solution 40. 

FIGURES 6 and 7 show another example of solutions obtained 
using embodiments of the coordinate transformation method 

20 and system of this invention. The model used to obtain the 
solutions shown in FIGURES 6 and 7 is a one-quarter segment 
of a well in the center of a circular reservoir model. The 
reservoir radius is 1,500 feet, and a relatively coarse mesh 
is used in the completion model, with an inner radius of .25 

25 feet, an outer radius of 100 feet, and having four 

intermediate computational nodes between these radii. As in 
the example of FIGURES 4 and 5, the reservoir is initially 
at 2,000 psi and a constant withdrawal is imposed at the 
well at all times (0-10 days) . 
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As can be seen in FIGURE 6, for radial distances beyond 
approximately 200 feet, the Cartesian operator solution 50, 
and the transform coordinate solution 60 track fairly 
closely to the analytical solution 40 for the pressure 
5 profile over the ten-day period- However, as shown in 
FIGURE 7, the near- well (here, out to 25 feet radial 
distance from the well bore) data clearly shows that the 
transform coordinate solution 60 is much more accurate in 
the near-well regions than the Cartesian operator solution 
10 50. FIGURE 7 is an expanded graph of the same solutions 

shown in FIGURE 6, for the near-well regions of .25 feet out 
to approximately 25 feet. In this near-well region, the 
transform coordinate solution 60 is a much better 
approximation of the observed analytical solution 40. 

15 Referring back to FIGURES 2 and 3, it can be seen that if 
any two points from FIGURE 2 are taken as one element, a 
large number of these points are necessary to create an 
approximation of the curve. By contrast, in the coordinate 
transformation space shown in FIGURE 3, the solution is 

20 linear and intermediate points are not necessary. Any two 
points can be used to pick up all variations and give exact 
solutions. 

The coordinate transformation of various embodiments of the 
present invention comprises a mathematical mapping from one 

25 coordinate space to another coordinate space and can be 
implemented as computer executable software instructions 
stored in memory within a computing device. The memory can 
be RAM or ROM, or a fixed media such as a hard drive or 
magnetic media as known to those in the art. The computing 

3 0 device can be any computer capable of executing the software 
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instructions, for example, a Windows UNIX or Linux based PC, 
or a Macintosh. The embodiments of the method and system of 
this invention can be described as a method for transforming 
coordinates to produce highly accurate simulations of 
5 reservoirs or other physical systems, in particular, near a 
point of interest such as a singularity (e.g., a well bore). 

Embodiments of the method and system of the present 
invention can be applied beyond the singularity of a well 
bore. Embodiments of this invention can be used to model 

10 and solve for properties of horizontal wells, fractures 
and/or faults. Embodiments of this invention are 
contemplated to be used in any situation where there is some 
potential for introducing a transformation .that more closely 
represents the behavior of an analytical solution. The 

15 transformation of the embodiments of the present invention 
is analogous to a case in solid mechanics for the modeling 
of heat transfer. 

Parameters other than pressure profiles can be modeled using 
the embodiments of the method and system of this invention. 

20 The embodiments of this invention comprise a method for 

modeling singularities of various types, with applications 
in many different industries. For example, point loads, 
temperature distributions, and heat sinks can be modeled 
using the embodiments of this invention. A reservoir could, 

25 for example, instead be imagined to be a metal plate, with a 
heat source applied, such as a flame or a laser beam, to 
heat it up. Mathematically, this situation is identical to 
modeling of a pressure profile in a reservoir. Any type of 
singularity can thus be modeled using embodiments of the 

30 method and system of this invention. 
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The embodiments of the method and system of this invention 
can provide faster, less computationally intense, and more 
accurate modeling and solutions than previous methods and 
systems for modeling physical systems near a singularity in 
5 that, unlike prior art methods and systems, they do not just 
apply standard finite element techniques in a Cartesian 
coordinate space. Instead, embodiments of the present 
invention apply coordinate transformation to the 
computational space. Transformation is applied and 
10 computation performed, and the results obtained are still in 
a Cartesian coordinate space. As far as the user is 
concerned, the method of the embodiments of the present 
invention is transparent. 

The embodiments of the method and system of this invention 
15 can be implemented within a user interface such as that 
disclosed in related U.S. Patent Application entitled 
"SYSTEM AND METHOD FOR DEFINING AND DISPLAYING A RESERVOIR 
MODEL" (the «Blitzer Application" ) . With the interface 
disclosed in the Blitzer application, for example, an 
20 oilfield engineer could specify the model parameters 

necessary to calculate the quantity that they are interested 
in, and the embodiments of this invention can perform the 
calculations in a manner transparent to the oilfield 
engineer . 

25 The embodiments of the method and system of this invention 
can be used to substantially overcome the inaccuracies of 
standard finite element analysis techniques when performing 
computations using a type of boundary condition representing 
the input or withdrawal of a physical quantity (for example, 

30 fluid or heat) from a singularity point. The embodiments of 
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the methods of this invention can comprise a method for 
performing a coordinate transformation in a computational 
space for more accurately calculating conditions around a 
singularity in a physical model, for example, of a reservoir 
5 or other physical system. 

Various embodiments of the method of this invention can be 
implemented in software and can be accessed through a user 
interface, such as the interface of the related Blitzer 
application, and can perform transformations in a method 
transparent to a user. These embodiments may include 
software or formal applications that embody methods such as 
those described herein, as well as computer readable media 
containing such applications and computers or other data 
processing systems that are configured to perform the 
described methods. Features of various embodiments of the 
present invention include transparent operation to the user, 
and the capability of allowing a user to set boundary 
conditions and variables of parameters that describe a 
system of interest. 

An advantage of the embodiments of the method of the present 
invention is that, unlike in traditional approaches for 
obtaining a finite element analysis solution near a 
singularity, fewer nodes are necessary to get the same or 
better accuracy. Previous methods and systems require 
specification of a fine spatial mesh around the singularity 
in order to obtain the same degree of accuracy possible with 
the embodiments of this invention. The embodiments of this 
invention thus permit using a coarser mesh without 
sacrificing accuracy over the prior art methods. Prior art 
methods, for example, might require 50 nodes to be specified 
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for the 6" to 20' radial distance out from a singularity to 
obtain a useful solution. Using the coordinate 
transformation method of the embodiments of this invention, 
nodes can instead be specified with much greater radial 
5 distance between them, for example, 10' to 20' between nodes 
(or even greater distances) . In fact, some embodiments of 
the coordinate transformation method of this invention can 
use nodes spaced 50' apart and still get as accurate a 
resolution as prior art methods where many more nodes must 
10 be specified between 6" and 50' out. 

The embodiments of the method and system of this invention 
work particularly well in the case of an isotropic system in 
which conductance in all directions is identical. An 
anisotropic system, where permeability (conductance) in the 
15 X and Y direction might be different, results in slightly 
less accurate solutions. Useable solutions are still 
possible in an anisotropic or heterogeneous system, but 
there may be extreme cases where the X direction 
permeability (conductance) maybe, for example, 1,000 times 
larger than the permeability in the Y direction, in which 
case the coordinate transformation solution of the 
embodiments of this invention can introduce a significant 



20 



error . 



Embodiments of the coordinate transformation method and 
25 system of this invention provide the ability to transform 
the equations necessary to solve a system to a form that 
looks Cartesian. Basically, the coordinate space is 
transformed such that it is very well represented by the 
finite element approximation. This is different from prior 
30 art techniques where an attempt is made to still use the 
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Cartesian equations, but to somehow embed a singularity into 
the finite element space. The embodiments of the method of 
this invention provide a technique that fixes the 
singularity (or removes the singularity) before the finite 
element method is ever applied to it, as opposed to trying 
to fix the finite element method by changing the way it 
operates. The underlying finite element analysis still 
exists on top of the transformation of the embodiments of 
the method of the present invention. 

Although the present invention has been described in detail 
herein with reference to the illustrative embodiments, it 
should be understood that the description is by way of 
example only and is not to be construed in a limiting sense. 
It is to be further understood, therefore, that numerous 
changes in the details of the embodiments of this invention 
and additional embodiments of this invention will be 
apparent to, and may be made by, persons of ordinary skill 
in the art having reference to this description. It is 
contemplated that all such changes and additional 
embodiments are within the spirit and true scope of this 
invention as claimed below. 
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T T Tg r*TiATMKD IS: 

1. A method for predicting the behavior of a physical 
system, comprising the steps of: 

creating an equation in a first coordinate system to model 
an aspect of said physical system; 

applying a coordinate transformation to said equation for 
transforming said equation from said first coordinate system 
into a second coordinate system more closely representative 
of an analytical solution to said equation; 
solving said equation in said second coordinate system to 
obtain a solution; and 

transforming said solution back to said first coordinate 
system. 

2. The method of Claim 1, further comprising the step of 
generating a mesh to model said physical system. 

3. The method of Claim 1, wherein said first coordinate 
system is a Cartesian coordinate system. 

4. The method of Claim 1, wherein said second coordinate 
system is a radial coordinate system. 

5. The method of Claim . 4, wherein said radial coordinate 
25 system comprises the coordinates (s, 6, z) , wherein s=ln(r) . 

6. The method of Claim 1, wherein said second coordinate 
system is a coordinate system where said solution is a 
linear solution. 
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7. The method of Claim 1, wherein said physical system is 
a well bore. 

8. The method of Claim 1, wherein said physical system is 
5 a reservoir having a singularity. 

9. The method of Claim 8, wherein said singularity is a 
well bore. 

10 10. The method of Claim 1, wherein said aspect is a 
physical property of said system. 

11. The method of Claim 10, wherein said physical property 
is pressure. 

15 

12. The method of Claim 1, wherein said physical system is 
a point load or a heatsink. 

13. The method of Claim 1, wherein said physical system 
20 comprises a fracture or fault. 

14. The method of Claim 1, further comprising the steps of: 
creating a second equation in said first coordinate system 
to model a second aspect of said physical system; 

25 solving said second equation in said first coordinate system 
to obtain a solution to said second equation; and 
combining said transformed solution to said first equation 
and said solution to said second equation in said first 
coordinate system to obtain a combined solution. 

30 



WO 02/03263 



27 



PCT/US01/20814 



15. The method of Claim 1, wherein said coordinate 
transformation comprises a mapping from said first 
coordinate system to said second coordinate system. 

5 16. The method of Claim 1, wherein said aspect is radial 
fluid flow near a singularity. 

17. The method of Claim 16, wherein said singularity is a 
well bore in a reservoir. 

10 

18. The method of Claim 1, wherein said coordinate 
transformation is performed via computer executable software 
instructions . 

15 19. The method of Claim 1, wherein said method is 
implemented using a graphical user interface. 

20. A computer- readable medium containing a plurality of 
instructions embodying a method for predicting the behavior 
20 of a physical system, said method comprising the steps of: 
creating an equation in a first coordinate system to model 
an aspect of said physical system; 

applying a coordinate transformation to said equation for 
transforming said equation from said first coordinate system 
25 into a second coordinate system more closely representative 
of an analytical solution to said equation; 
solving said equation in said second coordinate system to 
obtain a solution; and 

transforming said solution back to said first coordinate 
30 system. 
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21. The computer -readable medium of Claim 20, said method 
further comprising the step of generating a mesh to model 
said physical system. 

5 22. The computer- readable medium of Claim 20, wherein said 
first coordinate system is a Cartesian coordinate system. 

23. The computer-readable medium of Claim 20, wherein said 
second coordinate system is a radial coordinate system. 

10 

24. The computer-readable medium of Claim 23 , wherein said 
radial coordinate system comprises the coordinates (s f 9, 

z) , wherein s=ln(r) . 

15 25. The computer-readable medium of Claim 20, wherein said 
second coordinate system is a coordinate system where said 
solution is a linear solution. 

26. The computer-readable medium of Claim 20, wherein said 
20 physical system is a well bore* 

27. The computer-readable medium of Claim 20, wherein said 
physical system is a reservoir having a singularity. 

25 28. The computer-readable medium of Claim 21, wherein said 
singularity is a well bore. 



30 



29. The computer-readable medium of Claim 20, wherein said 
aspect is a physical property of said system. 
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30. The computer-readable medium of Claim 20, wherein said 
physical property is pressure. 

31. The computer-readable medium of Claim 20, wherein said 
5 physical system is a point load or a heat sink. 

32. The computer-readable medium of Claim 20, wherein said 
physical system comprises a fracture or fault . 

10 33. The computer -readable medium of Claim 20, wherein said 
method further comprises the steps of : 

creating a second equation in said first coordinate system 
to model a second aspect of said physical system; 
solving said second equation in said first coordinate system 
15 to obtain a solution to said second equation; and 

combining said transformed solution to said first equation 
and said solution to said second equation in said first 
coordinate system to obtain a combined solution. 

20 34. The computer- readable medium of Claim 20, wherein said 
coordinate transformation comprises a mapping from said 
first coordinate system to said second coordinate system. 

35. The computer- readable medium of Claim 20, wherein said 
25 aspect is radial fluid flow near a singularity . 

36. The computer- readable medium of Claim 35, wherein said 
singularity is a well bore in a reservoir. 

30 37. The computer-readable medium of Claim 20, wherein said 
method is implemented using a graphical user interface. 
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38. A method for predicting the behavior of a physical 
system, comprising the steps of: 

creating an equation in a first coordinate system to model 
an aspect of said physical system; 
5 mapping said equation from said first coordinate system into 
a second coordinate system more closely representative of an 
analytical solution to said equation; 

solving said equation in said second coordinate system to 

obtain a solution; 
10 mapping said solution back to said first coordinate system; 
creating a second equation in said first coordinate system 

to model a second aspect of said physical system; 

solving said second equation in said first coordinate system 

to obtain a solution to said second equation; and 
15 combining said mapped solution to said first equation and 

said solution to said second equation in said first 

coordinate system to obtain a combined solution. 

39. The method of Claim 38, further comprising the step of 
20 generating a mesh to model said physical system. 

40. The method of Claim 38, wherein said first coordinate 
system is a Cartesian coordinate system. 

25 41. The method of Claim 38, wherein said second coordinate 
system is a radial coordinate system. 

42. The method of Claim 41, wherein said radial coordinate 
system comprises the coordinates (s, 8, z), wherein s=ln(r) . 
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43. The method of Claim 38, wherein said second coordinate 
system is a coordinate system where said solution is a 
linear solution. 

44. The method of Claim 38, wherein said physical system is 
a well bore. 

45. The method of Claim 38, wherein said physical system is 
a reservoir having a singularity. 

46. The method of Claim 45, wherein said singularity is a 
well bore. 



WO 02/03263 



1/5 



PCT/US01/20814 




Fig. 1 



WO 02/03263 



2/5 



PCT/US01/20814 




WO 02/03263 



3/5 



PCTYUS01/20814 




WO 02/03263 PCT/USO 1/208 14 

4/5 



2000 



-I0OO 



nsoo 



JL170O 
2 1600 



13DO 



Well Brassure vs. Time 

•^O Time Interval* 











































































— 
















































1 1 1 1.] 
































**** 
















1*1 


I I 1 I 























a s e 



10 



WO 02/03263 



5/5 



PCT/USO 1/20814 




Pressure Profile at 10 days 




INTERNATIONAL SEARCH REPORT 



A CLASSIFICATION OF SUBJECT MATTER 

IPC 7 G01V11/00 



Intematlonal Application No 

PCT/US 01/20814 



According 10 International Patent Classification (IPC) or lo both n ational classification and IPC 
B. FIELDS SEARCHED 



Minimum documentation searched (classification system Wowed by Cassation symbols) 

IPC 7 G01V 



Documentation searched other than minimum 



"documentation to the extent that such documents are included in the fields searched 



■ Eleeronic data base consulted during the internation al search (name of data base and, where practieal, search terms used) 

EPO-Internal , WPI Data, INSPEC, COMPENDEX 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category" 



Citation of document, with indication, where appropriate, of the relevant passages 

OZDOGAN YILMAZ: "Seismic Data Processing" 
1991 SOCIETY OF EXPLORATION GEOPHYSICS , 
P.O. BOX 702740 TULSA, OK 74170-2740 
XP002189158 

page 514, paragraph 2 -page 515, paragraph 
1 

US 6 078 869 A (GUNASEKERA DAYAL L) 
20 June 2000 (2000-06-20) 
column 25, line 5 - line 65 
column 30, line 8 - line 28 

US 3 302 710 A (0DEH AZIZ S) 

7 February 1967 (1967-02-07) 

column 4, line 5 - line 20 

column 5, line 65 -column 6, line 70 

-/- 



Relevant lo claim No. 



1,20 



1,2,20, 
21,38,39 



1,20,38 



I Q(j F urtner documents are listed in the continuation of box C. 

[ • Special categories ot cited documents : 

■A' document defining the general state of the art which is not 

considered to be of particular relevance 
■E' earlier document but published on or after the international 

filing date 

V document which may throw doubts on priority dairnfc) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 
•O" document referring to an oral disclosure, use, exhibition or 
, other means 
'P' document published prior to the international filing date but 
later than the priority date claimed 

Date of the actual completion of the internalional search 

25 February 2002 

Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL - 2280 HV Rijswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 

Form PCT7ISA/210 (second sheet) (July 1992) 



Patent family members are listed in annex. 



■T- later document published after the internationa filing date 
or priority dale and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

■X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

■Y' document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 
Date of mailing of the international search report 



11/03/2002 



Authorized officer 



Lome, B 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 





Form PCT/ISA/210 {continuation of second sheet) (July 1992) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

Information on patent family members 



International Application No 

PCT/US 01/20814 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



US 6078869 



20-06-2000 



US 
CA 
FR 
GB 
NO 



6018497 A 
2228544 Al 
2760848 Al 
2322702 A ,B 
980821 A 



US 3302710 



07-02-1967 NONE 



25-01-2000 

27- 08-1998 
18-09-1998 
02-09-1998 

28- 08-1998 



Foim PCT/ISA/210 (patent family annex) (July 1992) 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
FADED TEXT OR DRAWING 

<E BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

-Q^COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

P?LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBITS) SUBMITTED ARE POOR QUALITY 

□ OTHER: 



IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



